Internal standard-based analysis of microarray data. Part 1: analysis of differential gene expressions

Genome-scale microarray experiments for comparative analysis of gene expressions produce massive amounts of information. Traditional statistical approaches fail to achieve the required accuracy in sensitivity and specificity of the analysis. Since the problem can be resolved neither by increasing the number of replicates nor by manipulating thresholds, one needs a novel approach to the analysis. This article describes methods to improve the power of microarray analyses by defining internal standards to characterize features of the biological system being studied and the technological processes underlying the microarray experiments. Applying these methods, internal standards are identified and then the obtained parameters are used to define (i) genes that are distinct in their expression from background; (ii) genes that are differentially expressed; and finally (iii) genes that have similar dynamical behavio

Statistical monitoring of weak spots for improvement of normalization and ratio estimates in microarrays

BACKGROUND: Several aspects of microarray data analysis are dependent on identification of genes expressed at or near the limits of detection. For example, regression-based normalization methods rely on the premise that most genes in compared samples are expressed at similar levels and therefore require accurate identification of nonexpressed genes (additive noise) so that they can be excluded from the normalization procedure. Moreover, key regulatory genes can maintain stringent control of a given response at low expression levels. If arbitrary cutoffs are used for distinguishing expressed from nonexpressed genes, some of these key regulatory genes may be unnecessarily excluded from the analysis. Unfortunately, no accurate method for differentiating additive noise from genes expressed at low levels is currently available. RESULTS: We developed a multistep procedure for analysis of mRNA expression data that robustly identifies the additive noise in a microarray experiment. This analysis is predicated on the fact that additive noise signals can be accurately identified by both distribution and statistical analysis. CONCLUSIONS: Identification of additive noise in this manner allows exclusion of noncorrelated weak signals from regression-based normalization of compared profiles thus maximizing the accuracy of these methods. Moreover, genes expressed at very low levels can be clearly identified due to the fact that their expression distribution is stable and distinguishable from the random pattern of additive noise

From microarray to biology: an integrated experimental, statistical and in silico analysis of how the extracellular matrix modulates the phenotype of cancer cells

A statistically robust and biologically-based approach for analysis of microarray data is described that integrates independent biological knowledge and data with a global F-test for finding genes of interest that minimizes the need for replicates when used for hypothesis generation. First, each microarray is normalized to its noise level around zero. The microarray dataset is then globally adjusted by robust linear regression. Second, genes of interest that capture significant responses to experimental conditions are selected by finding those that express significantly higher variance than those expressing only technical variability. Clustering expression data and identifying expression-independent properties of genes of interest including upstream transcriptional regulatory elements (TREs), ontologies and networks or pathways organizes the data into a biologically meaningful system. We demonstrate that when the number of genes of interest is inconveniently large, identifying a subset of "beacon genes" representing the largest changes will identify pathways or networks altered by biological manipulation. The entire dataset is then used to complete the picture outlined by the "beacon genes." This allow construction of a structured model of a system that can generate biologically testable hypotheses. We illustrate this approach by comparing cells cultured on plastic or an extracellular matrix which organizes a dataset of over 2,000 genes of interest from a genome wide scan of transcription. The resulting model was confirmed by comparing the predicted pattern of TREs with experimental determination of active transcription factors

A Comprehensive and Universal Method for Assessing the Performance of Differential Gene Expression Analyses

The number of methods for pre-processing and analysis of gene expression data continues to increase, often making it difficult to select the most appropriate approach. We present a simple procedure for comparative estimation of a variety of methods for microarray data pre-processing and analysis. Our approach is based on the use of real microarray data in which controlled fold changes are introduced into 20% of the data to provide a metric for comparison with the unmodified data. The data modifications can be easily applied to raw data measured with any technological platform and retains all the complex structures and statistical characteristics of the real-world data. The power of the method is illustrated by its application to the quantitative comparison of different methods of normalization and analysis of microarray data. Our results demonstrate that the method of controlled modifications of real experimental data provides a simple tool for assessing the performance of data preprocessing and analysis methods

Internal standard-based analysis of microarray data2—Analysis of functional associations between HVE-genes

In this work we apply the Internal Standard-based analytical approach that we described in an earlier communication and here we demonstrate experimental results on functional associations among the hypervariably-expressed genes (HVE-genes). Our working assumption was that those genetic components, which initiate the disease, involve HVE-genes for which the level of expression is undistinguishable among healthy individuals and individuals with pathology. We show that analysis of the functional associations of the HVE-genes is indeed suitable to revealing disease-specific differences. We show also that another possible exploit of HVE-genes for characterization of pathological alterations is by using multivariate classification methods. This in turn offers important clues on naturally occurring dynamic processes in the organism and is further used for dynamic discrimination of groups of compared samples. We conclude that our approach can uncover principally new collective differences that cannot be discerned by individual gene analysis

Molecular networks discriminating mouse bladder responses to intravesical bacillus Calmette-Guerin (BCG), LPS, and TNF-α

<p>Abstract</p> <p>Background</p> <p>Despite being a mainstay for treating superficial bladder carcinoma and a promising agent for interstitial cystitis, the precise mechanism of Bacillus Calmette-Guerin (BCG) remains poorly understood. It is particularly unclear whether BCG is capable of altering gene expression in the bladder target organ beyond its well-recognized pro-inflammatory effects and how this relates to its therapeutic efficacy. The objective of this study was to determine differentially expressed genes in the mouse bladder following chronic intravesical BCG therapy and to compare the results to non-specific pro inflammatory stimuli (LPS and TNF-α). For this purpose, C57BL/6 female mice received four weekly instillations of BCG, LPS, or TNF-α. Seven days after the last instillation, the urothelium along with the submucosa was removed from detrusor muscle and the RNA was extracted from both layers for cDNA array experiments. Microarray results were normalized by a robust regression analysis and only genes with an expression above a conditional threshold of 0.001 (3SD above background) were selected for analysis. Next, genes presenting a 3-fold ratio in regard to the control group were entered in Ingenuity Pathway Analysis (IPA) for a comparative analysis in order to determine genes specifically regulated by BCG, TNF-α, and LPS. In addition, the transcriptome was precipitated with an antibody against RNA polymerase II and real-time polymerase chain reaction assay (Q-PCR) was used to confirm some of the BCG-specific transcripts.</p> <p>Results</p> <p>Molecular networks of treatment-specific genes generated several hypotheses regarding the mode of action of BCG. BCG-specific genes involved small GTPases and BCG-specific networks overlapped with the following canonical signaling pathways: axonal guidance, B cell receptor, aryl hydrocarbon receptor, IL-6, PPAR, Wnt/β-catenin, and cAMP. In addition, a specific detrusor network expressed a high degree of overlap with the development of the lymphatic system. Interestingly, TNF-α-specific networks overlapped with the following canonical signaling pathways: PPAR, death receptor, and apoptosis. Finally, LPS-specific networks overlapped with the LPS/IL-1 mediated inhibition of RXR. Because NF-kappaB occupied a central position in several networks, we further determined whether this transcription factor was part of the responses to BCG. Electrophoretic mobility shift assays confirmed the participation of NF-kappaB in the mouse bladder responses to BCG. In addition, BCG treatment of a human urothelial cancer cell line (J82) also increased the binding activity of NF-kappaB, as determined by precipitation of the chromatin by a NF-kappaB-p65 antibody and Q-PCR of genes bearing a NF-kappaB consensus sequence. Next, we tested the hypothesis of whether small GTPases such as LRG-47 are involved in the uptake of BCG by the bladder urothelium.</p> <p>Conclusion</p> <p>As expected, BCG treatment induces the transcription of genes belonging to common pro-inflammatory networks. However, BCG also induces unique genes belonging to molecular networks involved in axonal guidance and lymphatic system development within the bladder target organ. In addition, NF-kappaB seems to play a predominant role in the bladder responses to BCG therapy. Finally, in intact urothelium, BCG-GFP internalizes in LRG-47-positive vesicles.</p> <p>These results provide a molecular framework for the further study of the involvement of immune and nervous systems in the bladder responses to BCG therapy.</p

A dynamic model of gene expression in monocytes reveals differences in immediate/early response genes between adult and neonatal cells

Neonatal monocytes display immaturity of numerous functions compared with adult cells. Gene expression arrays provide a promising tool for elucidating mechanisms underlying neonatal immune function. We used a well-established microarray to analyze differences between LPS-stimulated human cord blood and adult monocytes to create dynamic models for interactions to elucidate observed deficiencies in neonatal immune responses. We identified 168 genes that were differentially expressed between adult and cord monocytes after 45 min incubation with LPS. Of these genes, 95% (159 of 167) were over-expressed in adult relative to cord monocytes. Differentially expressed genes could be sorted into nine groups according to their kinetics of activation. Functional modelling suggested differences between adult and cord blood in the regulation of apoptosis, a finding confirmed using annexin binding assays. We conclude that kinetic studies of gene expression reveal potentially important differences in gene expression dynamics that may provide insight into neonatal innate immunity

Regulatory network of inflammation downstream of proteinase-activated receptors

BACKGROUND: Protease-activated receptors (PAR) are present in the urinary bladder, and their expression is altered in response to inflammation. PARs are a unique class of G protein-coupled that carry their own ligands, which remain cryptic until unmasked by proteolytic cleavage. Although the canonical signal transduction pathway downstream of PAR activation and coupling with various G proteins is known and leads to the rapid transcription of genes involved in inflammation, the effect of PAR activation on the downstream transcriptome is unknown. We have shown that intravesical administration of PAR-activating peptides leads to an inflammatory reaction characterized by edema and granulocyte infiltration. Moreover, the inflammatory response to intravesical instillation of known pro-inflammatory stimuli such as E. coli lipopolysaccharide (LPS), substance P (SP), and antigen was strongly attenuated by PAR1- and to a lesser extent by PAR2-deficiency. RESULTS: Here, cDNA array experiments determined inflammatory genes whose expression is dependent on PAR1 activation. For this purpose, we compared the alteration in gene expression in wild type and PAR1(-/- )mice induced by classical pro-inflammatory stimuli (LPS, SP, and antigen). 75 transcripts were considered to be dependent on PAR-1 activation and further annotated in silico by Ingenuity Pathways Analysis (IPA) and gene ontology (GO). Selected transcripts were target validated by quantitative PCR (Q-PCR). Among PAR1-dependent transcripts, the following have been implicated in the inflammatory process: b2m, ccl7, cd200, cd63, cdbpd, cfl1, dusp1, fkbp1a, fth1, hspb1, marcksl1, mmp2, myo5a, nfkbia, pax1, plaur, ppia, ptpn1, ptprcap, s100a10, sim2, and tnfaip2. However, a balanced response to signals of injury requires a transient cellular activation of a panel of genes together with inhibitory systems that temper the overwhelming inflammation. In this context, the activation of genes such as dusp1 and nfkbia seems to counter-balance the inflammatory response to PAR activation by limiting prolonged activation of p38 MAPK and increased cytokine production. In contrast, transcripts such as arf6 and dcnt1 that are involved in the mechanism of PAR re-sensitization would tend to perpetuate the inflammatory reaction in response to common pro-inflammatory stimuli. CONCLUSION: The combination of cDNA array results and genomic networks reveals an overriding participation of PAR1 in bladder inflammation, provides a working model for the involvement of downstream signaling, and evokes testable hypotheses regarding the transcriptome downstream of PAR1 activation. It remains to be determined whether or not mechanisms targeting PAR1 gene silencing or PAR1 blockade will ameliorate the clinical manifestation of cystitis

Trex1 regulates lysosomal biogenesis and interferon-independent activation of antiviral genes

Innate immune sensing of viral nucleic acids triggers type I interferon (IFN) production, which activates interferon-stimulated genes (ISGs) and directs a multifaceted antiviral response. ISGs can also be activated through IFN-independent pathways, although the precise mechanisms remain elusive. Here we found that the cytosolic exonuclease Trex1 regulates the activation of a subset of ISGs independently of IFN. Both Trex1−/− mouse and TREX1-mutant human cells express high levels of antiviral genes and are refractory to viral infections. The IFN-independent activation of antiviral genes in Trex1−/− cells requires STING, TBK1 and IRF3 and IRF7. We also found that Trex1-deficient cells display expanded lysosomal compartment, altered subcellular localization of the transcription factor EB (TFEB), and reduced mTORC1 activity. Together, our data identify Trex1 as a regulator of lysosomal biogenesis and IFN-independent activation of antiviral genes, and shows dysregulation of lysosomes can elicit innate immune responses